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Abstract 

In this paper we study a catalytically-activated A + A — ► reaction taking place 
on a one-dimensional regular lattice which is brought in contact with a reservoir of A 
particles. The A particles have a hard-core and undergo continuous exchanges with the 
reservoir, adsorbing onto the lattice or desorbing back to the reservoir. Some lattice 
sites possess special, catalytic properties, which induce an immediate reaction between 
two neighboring A particles as soon as at least one of them lands onto a catalytic site. 
We consider three situations for the spatial placement of the catalytic sites: regular, 
annealed random and quenched random. For all these cases we derive exact results for 
the partition function, and the disorder-averaged pressure per lattice site. We also present 
exact asymptotic results for the particles' mean density and the system's compressibility. 
The model studied here furnishes another example of a ID Ising-type system with random 
multisitc interactions which admits an exact solution. 
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1 This work is dedicated to the memory of our collaborator and colleague Professor Alexander A. Ovchin- 
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1 Introduction. 



Reactions involving particles which may recombine only when some third substance - the 
catalytic substrate - is present [1, 2], but otherwise remain chemically inactive, are ubiquitous 
in nature and also widely used in a variety of technological and industrial processes. Within 
the past two decades much effort has been put in understanding the peculiarities of such 
catalytically-activated reactions (CARs). In particular, considerable theoretical knowledge 
was gained from an extensive study of a particular reaction scheme - the CO-oxidation in the 
presence of metal surfaces with catalytic properties [3] (see also Ref.[4] for a recent review). 
Remarkably, Refs.[3] have substantiated the emergence of an essentially different behavior 
as compared to the predictions of the classical, formal-kinetics scheme and have shown that 
under certain conditions such collective phenomena as phase transitions or the formation of 
bifurcation patterns may take place [3]. Prior to these works on catalytic systems, anomalous 
behavior was amply demonstrated in other schemes [5-7], involving reactions on contact 
between two particles at any point of the reaction volume (i.e., "completely" catalytic sysems). 
It was realized [5-7] that the departure from the text-book, formal-kinetic predictions is due to 
many-particle effects, associated with fluctuations in the spatial distribution of the reacting 
species. This suggests that, similar to such "completely" catalytic reaction schemes, the 
behavior of the CARs may be influenced by many-particle effects. 

Apart from many-particle effects, the behavior of the CARs might be affected by the very 
structure of the catalytic substrate, which often cannot be considered as being a well-defined 
geometrical object, but represents rather an assembly of mobile or localized catalytic sites or 
islands, whose spatial distribution is complex [1]. Metallic catalysts, for instance, are often 
disordered compact aggregates, the building blocks of which are imperfect crystallites with 
broken faces, kinks and steps. Usually only the steps are active in promoting the reaction. 
In porous materials with convoluted surfaces, such as, e.g., silica, alumina or carbons, the 
effective catalytic substrate is also only a portion of the total surface area because of the 
selective participation of different sites in reaction. Finally, for liquid-phase CARs the catalyst 
can consist of active groups attached to polymer chains in solution. 

Such complex morphologies render the theoretical analysis difficult. There are only a 
few available studies which concern disordered substrates such as found in the CO-oxidation 
scheme; here the disorder is believed to affect mainly the particles' adsorption and desorption 
[8-16]. On the other hand, for the situations in which the spatial distribution of the catalyst 
is random, only empirical approaches have been used, based mostly on heuristic concepts 
of effective reaction order or on phenomenological generalizations of the formal-kinetic "law 
of mass action" (see, e.g., Rcfs.fl] and [2] for more details). The important outcome of 
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such descriptions is that they provide an evidence of the existing correlations between the 
morphology of the chemically reactive environment and reaction kinetic and steady-state 
properties. On the other hand, their shortcoming is that they do not explain the mechanisms 
underlying the anomalous kinetic and stationary behavior. In this regard, exact analytical 
solutions of even somewhat idealized or simplified models, are already highly desirable since 
such studies may provide an understanding of the effects of different factors on the properties 
of the CARs. 

In this paper we study the catalytically-activated annihilation A + A — ► reaction in a 
simple, one-dimensional model with different (regular or random) distributions of the cata- 
lysts, appropriate to the just mentioned situations with the catalytically-activated reactions 
assisted by the active groups attached to polymer chains. More specifically, we consider here 
the A + A — > reaction on a one-dimensional regular lattice which is brought in contact with 
a reservoir of A partilces. Some portion of the lattice sites (marked by crosses in Fig.l) pos- 
sesses special "catalytic" properties such that they induce an immediate reaction A + A — ► 0, 
when at least one of two neighboring adsorbed A particles sits on a catalytic site. In this case 
these two particles react and instantaneously leave the chain. 

In regard to the distribution of the catalytic sites, we focus here on three different situa- 
tions. First, we consider the case when the catalytic sites are placed periodically, forming a 
regular sublattice. Next, we turn to the disordered case. We analyse first the case of annealed 
disorder and furnish an exact solution. Lastly, we analyse the behavior in the most complex 
case of quenched disorder, for which situation an exact solution is also derived. 

We note finally that the kinetics of A + A — ► reactions involving diffusive A particles 
which react upon encounters on randomly placed catalytic sites has been discussed already 
in Refs.[20, 21] and [22], and a rather surprising behavior has been found, especially in low- 
dimensional systems. A much simplier equilibrium model of an A + A — ► reaction on a 
ID chain with randomly placed catalytic segments has been solved in Refs.[17] and [18], by 
noticing that here the average pressure per lattice site coincides exactly with the Lyapunov 
exponent of a product of random two-by-two matrices, obtained in Ref.[19]. Additionally, 
the steady-state properties of contact A + A — > reactions between diffusive A particles, or 
A + A — > reactions between immobile A particles reacting via long-range reaction proba- 
bilities in systems with external particles input have been presented in Refs.[23] and [24, 25], 
respectively, which analysis has revealed non-trivial ordering phenomena with anomalous in- 
put intensity dependence of the mean particle density. This anomalous behavior agrees with 
earlier experimental observations [26]. For completely catalytic ID systems, the kinetics of 
A + A — > reactions with immobile A particles undergoing cooperative desorption has been 
discussed in Refs.[27, 28] and [29]. Exact solutions for A + A — > reactions in ID completely 
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catalytic systems in which A particles perform conventional diffusive or subdiffusive motion 
have been presented in Refs.[30] and [31], respectively. 

The paper is structured as follows: In section 2 we formulate our model and introduce 
the basic notations. Next, in section 3 we focus on situations with regular placements of 
catalytic sites. We calculate exactly the partition function of the model, the pressure per 
site, and present as well explicit results for the particles' mean density in the thermodynamic 
limit from which we determine the compressibility of the system. In section 4 we study the 
behavior of the system in the case of random annealed distributions of the catalytic sites. 
We show that in this case the model reduces to a one-dimensional lattice gas with an ef- 
fective three-particle repulsive interaction. We develop a combinatorial formulation of the 
model which allows us to obtain an exact solution. We present thus an exact expression for 
the disorder-averaged pressure, as well as exact asymptotic expansions for the mean particle 
density and the compressibility. In section 5 we turn to the very complex situation where the 
random distribution of catalytic sites is quenched. For this case, averaging the logarithm of 
the partition function over the states of the quenched random variables describing the cat- 
alytic properties of lattice sites, we find that the problem reduces in finite lattices with a fixed 
number of catalytic sites to an exact enumeration of all possible interconnected clusters. The 
weights of such clusters are calculated exactly in terms of a certain combinatorial procedure, 
and we find eventually an exact expression for the disorder-averaged pressure in the quenched 
disorder case. Additionally, we evaluate exact asymptotic expansions for the mean-particle 
density and the compressibility. We show that the behavior of these properties differs substan- 
tially, depending on whether the disorder is annealed or quenched. In particular, in systems 
with annealed disorder the mean particle density tends to unity when the chemical potential 
jjL — > oo for any p < 1, where p is the mean density of catalytic sites. On the other hand, mean 
particle density in the quenched disorder case tends to a finite value (1— p + p 2 )/(l+p 2 ) < 1 
as /x — ► oo. As well, in the annealed disorder case the compressibility appears to be a non- 
monotonic function of p for any u, while in the quenched disorder case the compressibility 
shows a non-monotonic behavior dbS 8b function of p only for fjj ^ l^crit — 

0~ 1 \n(2), where 

is the temperature. Lastly, in section 6 we conclude with a summary of results. Details of 
intermediate calculations are summarized in the appendices A,B and C. 

2 The model 

Consider a one-dimensional, regular lattice containing N adsorption sites (Fig.l), which is 
brought in contact with a reservoir of identical, non-interacting, hard-core A particles - a vapor 
phase, maintained at a constant chemical potential \i. The A particles from the vapor phase 
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can adsorb onto vacant adsorption sites and desorb back to the reservoir. The occupation of 
the "i"-th adsorption site is described by the Boolean variable rij, such that 

i 1, if the "i"-th site is occupied, 
rii = < 

y 0, otherwise. 

For computational convenience, we also add two special, boundary sites i = and % = N + 1, 
and stipulate that these sites are always unoccupied, no = n^+i = 0. Suppose next that some 

• • VAPOR • ^ 

(a) 

i % i T T x % T x X X 1 

1 N N+1 

Figure 1: One-dimensional lattice of N adsorption sites in contact with a vapor phase. The filled 
circles denote A particles with hard-cores. The crosses denote the adsorption sites with catalytic 
properties, (a) denotes a "forbidden" particle configuration. The sites i = and i = N + 1 are always 
unoccupied and non-catalytic, i.e. no = n^+i = and (q = (n+i = 0. 



of the adsorption sites possess "catalytic" properties (crosses in Fig.l) in the sense that they 
induce an immediate reaction A + A — > between neighboring A particles; that is, if at least 
one of two neighboring adsorbed A particles sits on a catalytic site, these two particles react 
and instantaneously leave the chain (desorb back to the reservoir). To specify the catalytic 
sites, we introduce the quenched variable d, so that 







1, if the "i"-th site is catalytic, i = 1, 2, . . . , N, 
0, otherwise. 



The sites at the extremities of the chain are supposed to be non-catalytic, i.e. Co = and 
Cjv+i = 0- As for the segment [1, N], we will consider several possible ways of spatial placement 
of catalytic sites, namely, regular and random. 

For a given distribution of catalytic sites, the partition function ^at(C) of the system under 
study can be written as follows: 

N N 

Zn(C) = lim exp (fin y^rii- Xy^dm (n,_i + n i+1 )), (1) 

A— >oo — ' \ / 

{m} t=i t=i 
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where the summation extends over all possible configurations {n^}, fi is the chemical potential 
which accounts for the reservoir pressure and for the particles' preference for adsorption. The 
parameter A stands for the catalytic activity, which is here taken to be infinitely large. Such 
a choice implies that the reaction between two neighboring As, in which pair at least one of 
^4s sits on a catalytic site, takes place instantaneously. Note that here ^jv(C) is a functional 
of the configuration ( = {d}. 

We stop to note that an analysis of general reaction and diffusion problems was previously 
done using effective Hamiltonians of spin systems. Examples are, e.g., the works by Doi [32], 
Zeldovitch and Ovchinnikov [33], Alcaraz et al [34] and Simon [35]. Our procedure here is 
different in that we consider equilibrium systems and take the limit A — ► oo, which allows us 
to present the basic results in closed form. 

Now, taking into account that 



lim^exp ( - \Cini{m-i + ra i+ i)) = (l - Cz™i™i-i) (l ~ Ci™i™i+i)> ( 2 ) 

(V), (3) 



A— >oo 

and setting 

z = exp 



we can rewrite eq.(l) as: 



N 



Zn(() ^H^'f 1 - CWi-i) (l - C^+i) (4) 
{m} i=i 

Hence, any two neighboring sites i and i — 1 appear to be coupled by a factor (1 — n^ni-i) 
when at least one of these sites is catalytic. These coupling factors are depicted in Fig.l as 
arcs connecting neighboring sites and signify that configurations {rij} in which the occupation 
variables rii and rij_i assume simultaneously the value 1 are excluded. We introduce now the 
notion of "cluster", as being a set of sites, all connected to each other consecutively by arcs. 
Thus, a if-cluster contains K sites. Note that the boundary between adjacent clusters is 
given by a pair of two neighboring non-catalytic sites, i.e. by two consecutive variables Q 
and (i+i which are both equal to zero (see Fig.l). Now, the chain decomposes into disjunct 
clusters, and consequently, the partition function Zn(() factorizes into independent terms, 
such that each one depicts its corresponding cluster. 

It may be also instructive to rewrite eq.(l) in terms of the Ising-type spin variables <7j = 
(2rij — 1), such that Oi = ±1. In terms of these variables, Z]^(() of eq.(l) reads: 

N N N 

Zn(C) = A lim^ exp ([/3/j, - \p)— ) ^ exp ( ^ + ^ Ji^ai+^j (5) 

{(Ti} i=l »=1 
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where p denotes the mean density of catalytic sites, p = N 1 di while 

A*i = ^-^(2Ci + Ci-i + Ci+i), and J t = + Cm)- (6) 

Consequently, the model under study can be also thought of as a version of a one-dimensional 
Ising-type model with site-dependent magnetic fields and site-dependent couplings (see, e.g., 
some seminal works [19, 36, 37], as well Refs.[38] and [39] for a recent review). Note, however, 
that in our case both the fields and the couplings are non-local, and that the local energy 
£ij = —^i^i — JiCTidi+i, can assume several different values, depending on the occupation 
variables of the neighboring sites and on their catalytic properties. 

We close this section by mentioning the results of the conventional mean-field approach, 
which depicts the evolution of our system [1, 2]. Discarding correlations between the occupa- 
tion of neighboring sites, i.e. setting m = n, one writes the following balance equation 

n = -pJCn 2 + g ads (I - n) - g des n, (7) 

where the overdot denotes the time derivative, K, stands for the elementary reaction constant, 
p is the mean density of catalytic sites (introduced to account for the reduction in the reaction 
rate due to the partial chain coverage by catalytic sites) , while the second and the third terms 
on the rhs of eq.(7) correspond to the usual Langmuir adsorption/desorption events; here g ads 
and gdes are the adsorption and desorption rates, respectively, g a ds/gdes = z. 
Equation (7) has the following equilibrium solution 



{gads + gdes) I / ^pK-g, 




— \naas i naes) i ; ^t"^3ads , , i ,, 

n = Trip <\~, ; 72+ 1 - 1 r ( 8 ) 



2 ^ 1 V (gads + gdes 



In the limit /C — > oo, (which is the limit of interest in the present paper), and with p and g a d s 
kept finite one finds that n — ► 0; from eq.(8) n vanishes as n ~ g a ds/pfc- In the following 
we proceed to show that the actual behavior of the mean density of the A particles turns 
out to be very different from eq.(8); this is due to the emerging correlations between the A 
particles. Note also that in the limit K, = (suppressed reaction), one recovers from eq.(8) 
the classic Langmuir result n = z/(l + z) [1, 2]. 



3 Regular placement of catalytic sites. 



To fix the ideas, consider first a situation in which the catalytic sites are placed periodically, 
with period L, so that d obeys: 

vN — 1 



(i = 5(i,nL + 1), with n = 0, 1, . . . 



L 



(9) 
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where [x] denotes the integer part of the number x, and 5(k, m) is the Kroneker-delta symbol, 

dr I 1, if k = m, 



C T l+k 171 1 0, otherwise, 

where C stands for any closed contour which encircles the origin counterclockwise while 
(k,m) e 1?. 



VAPOR 



(a) 



1 N N+l 

Figure 2: Periodic placement of the catalytic sites with the period L — 4. (a) denotes a forbidden 
particle conhguration. 



We have now to distinguish between two situations: namely, when L > 3 and when L = 1 
or L = 2. In the former case, evidently, the factors (1 — Ci n i n i±i) m e Q-(4) are non-overlapping 
(see Fig.2); then the partition function decomposes into elementary three-clusters centered 
around each catalytic site and (possibly) into uncoupled, "free" sites, i.e. sites unaffected by 
any of the factors (1 — Ci n i n i±i)- On the other hand, in the cases L = 1 and L = 2 we deal 
with totally interconnected clusters, spanning the entire chain, as one can deduce from Fig. 3. 
In fact, the role of L = 1 and L = 2 is, chemically speaking, identical. 



• w VAPOR • ^ 

(a) (a) 

1 x T x T x T x 1 x 1 — ' 

1 N N+l 

Figure 3: Periodic placement of the catalytic sites with the period L = 2. In this case all sites are 
coupled by factors (1 — nj7ij_i) and hence the occupation variables of any two neighboring sites can 
not assume the value 1 simultaneously; (a) denotes such "forbidden" particle conhgurations. 
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3.1 Periodic case with L > 3. 

In the case L > 3 the partition function in eq.(4) decomposes into the product 

Z N (() = 4 eff) (L) = Z«> zf\ 



(10) 



where the superscript "(reg)" signifies that we deal with the regular case, and where Zr, 
K = 1, K = 2 and K = 3, are the partition functions of one-, two- and three-clusters, while 
Mr, K = 1,2, 3, stands for the numbers of such clusters in the TV-chain, respectively. Now: 



IN - 1- 


-'(V. 


-N - In 


l L - 




- L J 



); A/" 2 = l + <5( 



TV- 1 



TV - 1 



L 



and 



//V- 1 


rJV- 1- 


)-* 


\N- 1] 


V L ' 


L L \ 


- L \ 



Ni = N -2 + 5( 
where we have used the evident "conservation" law 

3AA + 2AA2+M = 7V. 



(11) 
(12) 

(13) 



Note, however, that the number of the two-clusters is not extensive, i.e., it does not grow 
with TV: Such clusters can be present only on the boundaries, i.e., A/2 = 2 when (N — 1)/L is 
an integer and A/2 = 1 otherwise. Moreover, 

Z i= E ^ = (l + z); Z 2 = z ni+n2 (l-nm 2 ) = (l + 2z), (14) 

{ni=0,l} {ni,n 2 =0,l} 

and 

Z 3 = z ni+n2+n3 (l-n 1 n 2 )(l-n 2 n 3 ) = (l+3z + z 2 ). (15) 

{ni,n2,n3=0,l} 

Consequently, in the case of a regular, periodic placement of catalytic sites with period L, 
L > 3, we can calculate the pressure p( re9 \L) per site from the relation: 



(3P i - re9) {L) = lim 

N^oo 



lnZ^\L) 
N 



(16) 



which yields: 



/3P (re9 \L) 



/A/3 A/2 A/i 

lim — In Z 3 - — In Z 2 - — In Z\ 

N^oo \ N N N 

p ln(l + 3z + z 2 ) + (1 - 3p) ln(l + z), 



(17) 



where p = 1/L is the density of catalytic sites. 

The averaged density of adsorbed particles obeys 



n {re9 \L) 



1 

N 



N 



n, 



(18) 
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and hence, in the limit N — > oo it follows that 

„^) (i) = (1 _ 3 p) r i- + p I |±^. (,., 



We stop to note that in the limit p = (i.e. that of a completely non-catalytic chain) eq.(19) 
reduces to the standard Langmuir result. On the other hand, when p = 1/3 the number of 
free-sites in the periodic chain vanishes and the Langmuir contribution gets equal to zero; in 
this case the chain is composed totally of three-clusters. In this case, for very large z the 
occupation of lattice sites by adsorbed molecules equals 2/3, as it should. 
Finally, we have that the compressibility kx, defined as 

(3^ 1 k ^ 9n z dn 
n 2 dfi n 2 dz 

obeys 

rt-i t ("*>m - 1 ± 6z + llz2 + fc3 + zA ~ pz{ * + 8z + 19z2) (ou 

P kT {L) ~ z(l + 3z + Z 2 -4 P Z- P Z 2 ) 2 ■ W 

For p = the last equation reduces to the standard Langmuir result, k T an ^ = j3/z. 
3.2 Periodic cases with L = 1 or L = 2. 

We note first that since (1 — njrij+i) 2 = (1 — njTij+i), we evidently have that 



N 



Z N ee Z%^{L = 1) = Z<™\l = 2)= J2 zZr=in ' Hi 1 - 

{ni},i=l,...,N i=l 



• (22) 



n N+1 =0 

We hence conclude that the L = 1 and the L = 2 periodic systems are chemically equivalent. 

Now, we focus on the derivation of the form of Z^ according to eq.(22). To do this, we 
proceed as follows: performing first the averaging over the states of the variable n^, we have 
that Zjy can be written as 

N-2 

z N = z N _ x + Z yi zE ^ 1 th n i 1 - n ^+o i 1 - > ( 23 ) 

{m},i=l,...,N-l i=l 

where Zn-i obeys 

JV-1 



Z N -i = Yl zZZi 1 ** 11 (l - n t n i+l ) 
{m},i=l,...,N-l i=l 



(24) 

n N =0 

and is hence the partition function of a chain consisting of N — 1 sites. We notice next that 
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the second term on the rhs of eq.(23) vanishes when un-i = 1; thus it can be written as 

N-2 

{ni},i=i,...,JV-i i=i 

AT-2 

^ JJ (l_ nin . +1 ) , (25) 

{n i },i=l,...,JV-2 i=l njv _ i=Q 

where the expression on the right-hand-side is, evidently, the partition function of a chain 
containing N — 2 sites. Consequently, the partition function of an JV-site chain obeys the 
following recursion: 

Z N = Zat_i + zZ N - 2 (26) 

whose first three terms are given by eqs.(14) and (15). 

Next, to determine explicitly for arbitrary N we introduce the following generating 
function 

oo 

H(r) = J2 Znt n . (27) 

N=l 

Then, multiplying both sides of eq.(26) by t n ~ 2 and performing the summation, we obtain 
the following explicit result for H(t): 

H(t) = [zT 2 + (1 + Z)r)(l - T - = 

(zt 2 + (1 + z)t) 1 



■(-^+— ), m 

\ T + Tl To — T / 



y/1 + 4z V r + ri r 2 - r / 
where t\ and r 2 are given by 

n = Y z (ViT^ + l) (29) 

and 

T2 = ^(v / T+4^-l). (30) 

Next, expanding the terms in the second line on the rhs of eq.(28) in a Taylor series in powers 
of r and gathering terms entering with the same power, we find that Z^ obeys 

ZN - 2Vl + 4*r 2 * £jV ' (31) 

where 

V v ; (i + 2« + vT+4i) Wi/ / v ; 

We recall that the Z/v are, of course, polynomial functions of the activity z; expanding the 
rhs of eq.(31) in powers of z, we get 

[^ +1 )/ 2 1 /AT / + 1\ 



Z=0 
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where (^) denote the binomial coefficients, 

N\ _ j N\/U(N - Z)!,for N > I, 
I ) y 0, otherwise. 

For large values of z, it might be more convenient to use another representation; subtracting 
in eq.(33) z^ N+l ^ 2 and summing up the remaining terms, we find that Zn are given by 

Z N = Z ( N+i y 2 F N+2 (l/^), (34) 

where F n (x) are the Fibonacci polynomials [40], defined explicitly by 

K»»-i)/2] , _ , _ 1X 
Fn(x)= £ J (35) 

Finally, noticing that for z < oo one has (t-z/tx) < 1 and hence, that (t2/ti) n vanishes 
exponentially with N as N — > oo, we find in the cases when L = 1 and L = 2 that the pressure 
per site obeys, using eq.(16): 

(3p(^9)(L = lor 2) = ln(^- 1+ ^ + 1 ), (36) 
while the average density in an infinite chain is, using eq.(18): 

n (reg) = 1 Q r 2) = 1 2 ^ - . (37) 

1 + Az - VI + 4z 

In the limit z = oo, the roots 7~i and T2 are equal, Ti = T2, this signals the emergence of long- 
range order. Actually, in this case n^ re9 ^(L = lor 2) = 1/2 and the particles' distribution on 
the lattice is periodic. 

Finally, to close this section, we derive the compressibility of the particle phase. In the 
thermodynamic limit, for regular placement of the catalytic sites with the period L = 1 or 
L = 2 it obeys: 

p-'k^iL = lor 2) = ; -. (38) 

T V ' v / TT4l(l + 2z - 

Note that contrary to the expression in eq.(21) which holds for L > 3, here k^ e9 \L = lor 2) 
is a non-analytic function of the activity z when z — > 00. 

4 Random placement of catalytic sites. Annealed disorder. 

Consider next the situation with annealed disorder, which is realized, for instance, when the 
catalytic property may move (diffuse) very quickly. In this case, the disorder-average pressure 
p( ann \p) per site and the average density are given by, respectively, 

P (-) (p)= Hm ^P-, (39) 
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and 

' N 



n (ann) {p) = n . = ^^p(«nn )(p)j (4Q) 

i=l 



where Zjv(C) is again the partition function of eq.(l). 

We can now average directly the partition function in eq.(l) over the placement of the 
sites with catalytic property: 

N I N \ 

{Z N (0) = lirn^ ^ exp ^ ^ / exp ( - A ^ Ci™i(™i-i + ) = 

"^{n;} i=l \ i=l / 

iV V 

= lim^ ^ exp (/fy ^ J| ^exp ( - XQn^rii-i + n i+ i)) ^ = 

{«,} »=1 i=l 

= ^z^i^lim^ Yl (pexp ( - An^-i +n i+1 fj + 1 -p). (41) 

{n;} i=l 

Noticing next that 

,. / ( x , , ^ , r \ / if tti(ni_i +n m ) > 0, 

hm pexp -Arijni-i+fii+i +l-p =< 

A— »oo ; ; | l,ifn i (n i _i+n i+1 ) = 0, 

and hence, that 

lim (pexp(-An i (n i _i+n m )) + = (1 - p) n ^ 1 ~ ^ ~ "i-iX 1 ~ (42) 

we find that the averaged partition function in eq.(l) attains the form 

(^(0) = Y J ^ ni ^~P) Lt ^\ (43) 



where <I>j is a Boolean variable: 

= (riirii+i + riim-i - ni-inin i+ i 



I 1, if nj(nj_i + n i+ i) > 0, 
[ 0, if rn{ni-i + n i+ i) = 0. 

Evidently, ^ is non-local and depends on the environment of the "i"-th site. As a matter of 
fact, the local energy q at the "i"-th site, q = — fliirii — ln(l— p)(njnj_|_i+njnj_i — rij_injnj + i), 
assumes 8 different values depending on nj_i, rij and rtj+i. Before we proceed further, it might 
be also instructive to rewrite eq.(43) in terms of the spin variables crj. Then eq.(43) takes the 
form 

3 N N N 

(Z N (()) = exp((/3^ + -ln(l-p))-)^exp(/3 A i / ^ C j J + Ji^^ +1 - 

{o-j} i=l i=l 

N N 

- J2 ^ Ci-iOt+i - J3 ^ Q'i-iO'iO'i+i) , (44) 
i=i i=i 
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where 

M' = /i + ^In(l-P); Ji = ^ r ^<0; and J 2 = J 3 = < 0. (45) 

As one may notice, the expression in eq.(44) represents a combination of two well-known 
Ising-type models: namely, the first three terms in the exponent define the antiferromag- 
netic ANNNI (axial next-nearest neighbor Ising) chain [41]. In our case, the competition 
ratio k = J2/J1 equals k = 1/2, which is, in fact, a non-trivial special point (the so-called 
multiphase point [42]) of the ANNNI model. On the other hand, the fourth term in the 
exponent corresponds to the three-spin interaction model [43]. Both models have been ex- 
tensively studied within the last two decades in various contexts [42, 43] and show interesting 
equilibrium and dynamic properties, see e.g., Refs.[44, 45]. We are, however, unaware of an 
exact solution of the one-dimensional combined model of eq.(44). Below we will furnish such 
an exact solution using a combinatorial approach. 

Note now that *j always equals zero for unoccupied sites, (rii = 0), but attains the 
value = 1 only for occupied sites, (rii = 1), which have at least one (or two) occupied 
neighboring sites (see Fig.4). Otherwise, for isolated occupied sites (elementary sequences 



o 

I 



• o 

J l_ 



• o • o 

J I I l_ 



• o 

J l_ 



• OOO 

J l l l_ 



Figure 4: Values of the variable corresponding to a given configuration {rii}. 



with (rii-i = 0,rij = l,nj+i = 0)) the variable ^ equals 0. Consequently, for any given 
realization {rii}, one has that 

N 

J2^i=^+[{ni}]-Ks[{ni}}, (46) 
i=i 

where A/+[{nj}] is the number of lattice sites on which (in a given realization {rii}) the occu- 
pation variable rii assumes the value 1, while A/i s [{nj}] is the realization-dependent number 
of isolated occupied sites (elementary cells of the form (0, 1, 0)). Hence, the partition function 
in eq.(43) can be rewritten as 

(Zn(0) = £ (,(1 -p)f +[K}1 (l -p)-M.M (47) 

{rii} 
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Next, ordering the entire set of 2^ different realizations {rij} with respect to the total number 
of adsorbed particles which they contain, i.e. A/+[{rij}], we can recast eq.(47) into the form: 



N n+ N - N++1 



(Zn(Q) = £ (z(l-p)) + £ (l-Pr m M m (N + \N), (48) 

N + =0 m=0 

where M m (N + \N) stands for the number of realizations {rij} that have a fixed N + and contain 
exactly m elementary cells (0, 1,0). 

4.1 Calculation of M m (N + \N). 

To evaluate M m (N + \N) we now proceed as follows. Let us consider a given realization {n^} 
with N + filled (and 7V_ = N — N + vacant) sites and specify the lattice positions of the vacant 
sites by introducing a set of intervals {lj}, j = 1, . . . , iV_ + 1, such that the interval l\ connects 
a boundary site i = with the first vacant site, the interval h connects the first vacant site 
with the second one, and etc, while the last interval In-+i connects the last vacant site with 
the boundary site i = N + 1 (see Fig. 5). 

These intervals, which uniquely define the positions of the vacant sites in each given 
realization {rij}, obey the "conservation law": 

h + h + k + • • • + In.+i = N + 1 (49) 

Then, since each interval containing exactly two lattice units corresponds to an isolated 
occupied site, M m (N + \N) equals the number of different solutions of eq.(49), constrained by 
the condition that in each sequence {k} obeying eq.(49) m of the iV_ + 1 intervals are equal 
to 2, while the rest can assume any value except 2. Hence, M m (N + \N) can be represented as 

M m {N + \N) = ( N ~ + ^ P m (N + \N). (50) 

Here the binomial coefficient accounts for all possible choices of m intervals from the set of 
iV_ + 1 intervals, while P m (N + \N) stands for the number of different solutions of the equation 

h + h + k + ■■■ + lN.-m+i = N-2m + l, (51) 

where each of intervals 1^, k = 1,2, ... , iV_ — m + 1, can assume one of the values 4 = 
1, 3, 4, 5, — Making use of the integral representation of the lattice delta-function in eq.(10), 
we can write P m (N + \N) as 



(52) 
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f t t 



_L 



1 



N N+l 



Figure 5: Definition of a given configuration {ni} in terms of the intervals {4} connecting sequentially 
unoccupied lattice sites. 



where the hat above the summation symbol signifies that summation runs over all possible 
values Ik = 1, 3, . . . , oo, excluding the value Ik = 2. Performing the summation, we find 



(jV+ - rnj 



N-N + +l-m 



where 



S r =((l-r)- 1 -r). 



(53) 
(54) 



4.2 Explicit form of (Z N (()). 

Substituting eqs.(53) and (50) into eq.(48), we find that the partition function obeys 

N 



JL N+=0 T m=0 v 



which yields, upon summing over m and N + , the following result 
2(1 —p)( 1 _ 3V ~Q cos arccos (r/ \/-<5 3 



(Zn(C)) 



+ 



z(l-p) 



(55) 



9 (~Q)P X l + 2cos (| arccos (fl/v^Q 5 )) V ' ] 
1 + 3 v /3 Qsin (-arcsin (i?/\/-Q 3 )) / Z (i_ p )^+2 



Af+2 



+ 



2 sin Q + - arccos (Vv 73 ^)) - 1 
1 + 3 v /Z Qsin ^arcsin (-R/V-Q 3 )) - p)\ JV+a-j 
1 - 2 sin - | arccos (r/ y/=Q*) y h ' 



where R and Q are auxiliary functions, which obey 

I _ (1-p) 
9 

and 

1 , l/(l-p) 



3p 



(l + z(l-p)) < 0, 



27 + 6 



(l + z(l-p))- 



3z(l -p) 2 



(56) 



(57) 



(58) 



p p 
Equation (56) defines (Zjv(C)) for arbitrary values of p, z and chain's length N. The derivation 
of eq.(56) is presented in Appendix A. 
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4.3 The thermodynamic limit N — > oo. 

We turn next to the thermodynamic limit aiming to calculate the disorder-averaged pressure 
per site in the annealed disorder case and the mean density of adsorbed particles. In the 
limit N — > oo only the smallest root in absolute value matters. To select the appropriate 
root, it suffices to plot the combination (pi — l/S)/2y/—Q, i = 1,2,3, versus the variable 
x = Rj \/—Q 3 , which is defined on the interval [—1,1]. This plot is presented in Fig. 6 and 
shows that the smallest root is t = ti. 




Figure 6: Plot of (U - l/3)/2y/=Q, i = 1,2,3, versus the variable X = R/^Z-Q 3 . The dotted line 
(1) gives (ii - l/3)/2V^Q, the solid line (2) gives (i 2 - l/S)/2y/=Q, and the dashed line (3) gives 
(i3 — l/3)/2\/— Q, respectively. 



Consequently, taking into account that 

iresin (r/^-Q 3 )) 

> 0, (59) 



1 + 3-v 7 — Q sin ^- arcsin (r/ \/—Q' 



2 sin + - arccos (i?/\/-Q 3 )) 
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is positive definite, we find from eq.(39) that the pressure per site is given by 

/ (l - 6V=Qsin arcsin (ll/y/^Q^^ ^ 



In 



3z(l-p) 



(60) 



which expression determines p( ann )(p) for arbitrary values of the parameters p and 2. 

Consider now the form of p( aran )(p) in the limits p — > and p ~ 1. In the limit of a 
vanishingly small concentration of catalytic sites, we have from eq.(60) that: 

ppiann) {p) = ln(1 + z) _ j^^p + (p 2 ) . (61) 

Note that the first term in the last equation is, as it should be, just the standard Langmuir 
adsorption isotherm. 

In the case p ~ 1, we expand first t 2 /z(l — p) in powers of (1 — p). This yields, 
t 2 vT+4z-l (1 + ^)^1 + 4^-32-1 



■(l-p) 2 + 0((l-p) 3 ), (62) 



z(l-p) 2z 2zVl + 4z 

which expansion, as one can check by comparing the first and the second terms in eq.(62), 
makes sense only when (1 — p) <C (4/2) 1 / 4 . As a matter of fact, p = 1 has here a special role, 
as will be shown in the following. Then, we find that the pressure per site is 

0P t«nn) {p) = 1q (VT+Tz + U 4 2 2 + 5z - -^JY+Tz + 1 - VT+Tz ({l-pf) 

y ' v 2 ) (yr^4j_ i)(i + 4z ) v F > v 11 ) 

(63) 

Note now that the first term in eq.(63), which defines the disorder-averaged pressure per 
site in the limit p = 1, coincides exactly with the result we obtained earlier in the case of a 
regular placement of catalytic sites with the period L = 1 (or L = 2), eq.(36). As a matter 
of fact, one could notice that the partition function in the annealed-disorder case, eq.(48), 
will coincide for p = 1 with the partition function for the regular case, eqs.(31) and (32), by 
just analysing the behavior of the coefficients M m (N + \N), eqs.(50) and (53). To show this, 
we note first that from eqs.(50) and (53) one has that M m (A + |A^) = for m > N + . This is, 
of course, quite evident, if we recall that M m (N+\N) stands for the number of realizations 
{rii} having a fixed number N + of adsorbed particles and a fixed number m of elementary 
cells containing one adsorbed particle surrounded by two vacant sites: hence, the number 
of realizations {n{\ in which m exceeds N + equals zero. Further on, one notices that when 
p = 1, in eq.(48) only the terms M m=N+ (N+\N) with N + < N_ + 1 = N + 1 - N + matter. 
Hence eq.(48) becomes 

[(JV+l)/2] 

(Z N (Q)= E M N+ (N + \N)z N \ (64) 

N + =0 
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where 



Ic r - 



2iri \ N + J Jc t 



1 [N_ + 1 - N+\ f dr / ^ v N-+1-N+ (N_ + 1 - N. 



« ) s ( W+ ) < 65 > 



On comparing eqs.(64) and (65) with eq.(33), we have in the limit p = 1 that the partition 
function, (and hence, the pressure per site), in the annealed disorder case coincides with the 
partition function of a chain on which the catalytic sites are placed regularly with period 
L = 1. 

Finally, differentiating eq.(60) with respect to the variable z and making use of eq.(40), 
we find that in the annealed disorder case the average density of adsorbed particles obeys: 

n (ann) {v) = j + 3z(l - p) 2 ^ 

1 - Gv^Qsin (- arcsin (i?/V-Q 3 )) 
Asm (1 arcsin (r/ ^Q*)) + + ^ 2/g3 cos arcsin (^/v 73 ^))] , (66) 



where 



and 



A = (p(3zp 2 - 2p(l + 3z) + 3(1 + z))) , (67) 



B= 3 / (2V+p(5-4 Z ) + 2z-7)^ j (6g) 



(3zp 2 - 2p(l + 3z) + 3(1 + z)) 

Note that the result in eq.(66) differs from the mean-field prediction n = 0, which follows 
from eq.(8) for instantaneous reactions, /C = oo. 

In Fig.8 we present a plot of n^ ann ^ (p) versus p for different values of the activity z. In Fig.8 
we also compare the behavior of n^ ann \p) with the behavior found in the case of quenched 
disorder (see the next section). In the limits p <C 1 and p ~ 1, we find that n^ ann \p) follows, 
respectively, 

and 

W = 1 - 1 + 4 ,-,/raz + (TT^ (1 - p)2 + ((1 - p)3) ' (70) 

Note that the first term in eq.(69) is, as it should be, just the Langmuir adsorbtion isotherm, 
while the first term in eq.(70) coincides with our earlier result for the average density of 
adsorbed particles on the periodic chain with L = 1 or L = 2, see eq.(37). 
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Lastly, we analyse the behavior of the compressibility nn ^ (p) for the annealed disorder 
situation: 

g-^nn) (v) _ z dn^\p) 
We find that nn ^ (p) shows the following asymptotical behavior: For p <C 1 we have 



P- 1 4 nn \p) = l + f^P + 0( P >), (72) 



while for p ~ 1 the compressibility obeys 

0~^ nn \p) = ^kp 9 \L = lor 2) + (1+ 4 ^ )3/2 (l-p) 2 + 0((1 - P) 3 ). (73) 

Here k^ e9 \L = lor 2) is as previously defined, eq.(38), and represents the compressibility of 
a completely catalytic chain. Note that the result in eq.(73) signifies that in the annealed 
disorder case the compressibility is a non-monotonic function of the mean density p of catalytic 
sites. This can be seen immediately if one notices that, first, for any fixed z, one has kp 9) (L = 
lor 2) > /cy Cm - ) = f3/z, (or in other words, that for any fixed z the compressibility of a 
non-catalytic (Langmuir) system is always smaller than the compressibility of a completely 
catalytic system) and second, that k^ nn ^ (p) approaches kf^ 6£ ^ (L = 1 or 2) from above, since 
the function 4z 2 /(l + 4z) 3 / 2 is always positive. 

We close this section with some comments concerning the large- z behavior of p( ann \p) 
and n ( - ann \p). As a matter of fact, it appears that in the annealed disorder case the large- 
z behavior of p( ann \p) and n t - ann \p) for p arbitrarily close but strictly less than unity is 
completely different from the large- z behavior of these parameters in the case when p = 1. 
This implies, of course, that p = 1 is a special point. More specifically, we find that for 
z » (1 — p)~ 2 the disorder-averaged pressure per site obeys 

PP^HP) = laM + Ml - P) + ^ - ii^g, + O(i) , (74) 
which implies that the mean density follows: 

» ( - , W = 1 -7rh5 + OT + °(7)- (75 » 



This asymptotic behavior should be contrasted to the asymptotic behavior which holds in the 
p = 1 case. When z — > oo, we find from eqs.(36) and (37) the following results: 

pp(reg) ( L = 1 or 2) = - ln(z) + — L_ 3— + - + f J-") , (76) 

V ' 2 y ' 2zV2 48 z 3 / 2 1280 z 5 / 2 \z 7 / 2 ) V ; 



and 



n ( re s)(L = lor 2) = 1 TT r + — ^^ + 0[-^\. (77) 

V ; 2 4zV2 32 z 3 / 2 512 z 5 / 2 V ^ 7/2 
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This signifies, in particular, that for p arbitrarily close but not equal to unity, the mean 
density is equal to 1 as z = oo, while for p strictly equal to unity the mean density n (re9 \L = 
lor 2) = 1/2. On physical grounds, such a behavior can be understood as follows. In the 
annealed disorder case, instead of averaging the logarithm of the partition function in eq.(l), 
we average the partition function itself and thus operate with an effective, "annealed" partition 
function in eq. (43) . Here, the strict constraint that no two particles can occupy simultaneously 
neighboring sites if at least one of them sits on the catalytic site, is replaced by a more tolerant 
condition (see, eq.(44)), which allows for such pairs to be present at any site, while a penalty 
of 2 ln(l — p) is to be paid. For any finite p < 1 such a penalty can be overpassed by increasing 
the chemical potential. Thus for j3fi » — 21n(l — p) (or, equivalently, for z>(l- p)~ 2 ) one 
expects essentially the same behavior regardless of the value of p, and, in particular, that 
n( ann \p) — > 1 as z — ► oo. On the other hand, for p = 1 the penalty for having a pair of 
particles occupying neighboring sites becomes infinitely large and can not be compensated 
by any increase of the chemical potential. The behavior of n^ ann \p) as a function of z for 
different values of p is depicted in Fig.9 and compared with the behavior obtained for it in 
the case of quenched disorder. 



5 Random placement of catalytic sites. Quenched disorder. 

We finally turn to the most challenging situation - the case of quenched randomness in the 
placement of catalytically active sites. We begin by introducing one auxiliary function. Con- 
sider a chain of length N which contains a fixed number ./V — N nc of catalytic and hence, 
a fixed number N nc of non-catalytic sites, the latter being placed at the positions {X n }, 
n = 1, 2, . . . , N nc - We denote the partition function of such a chain as Z^({X n }). Evidently, 
Zn({X u }) obeys eq.(l) and eq.(4) with 

f 0,ifthei£{X n }, 
[ 1, otherwise. 

Then, the logarithm of the partition function in eq.(l), averaged over all realizations of the 
quenched random variable {Q}, can be formally written as 

N 

(lnZ N (()) = P N ~ Nnc (l-p) Nnc £m^({*n}), (78) 

N nc =0 {X n } 

where the sum with the subscript {X n } signifies that the summation extends over all possible 
placements of N nc non-catalytic sites on an iV-chain. 

Next, similarly to the approach used in the previous section, we introduce a set of N nc + 1 
intervals {l n } determined by consecutive non-catalytic sites, such that l n = X n — X n _i (with 
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Xq = 0) and In„ c +i = N + I — X^ nc . That is, the first interval extends from the boundary 
(non-catalytic, unoccupied) site i = to its closest non-catalytic neighboring site, the second 
interval extends from this non-catalytic site to the following one, and so forth, while the 
interval /Ar nc +i goes from the last non-catalytic site inside the chain to the boundary site 
i = N + 1. In terms of these intervals, eq.(78) can be rewritten as 

N 

(\nz N (o) = P N - Nnc (i-p) Nnc J2 lnZ "tt 1 ^ ( 79 ) 

N nc =0 {/„} 

where Zj^({l n }) stands now for the partition function Z]y({X n }) in which we have just ex- 
pressed the positions of the non-catalytic sites using the set {l n }, while the summation with 
the sign {l n } denotes now the summation over all possible solutions of the equation, analogous 
to eq.(49), 

h + h + h + ■ ■ ■ + l Nnc +i =N+l, (80) 

where each k> 1. 

For each given set {l n } of intervals we have that the partition function of an iV-chain 
decomposes into that of smaller clusters, 

Z N ({l n }) = zf li{ln}lN) z^ 2i{ln}lN) zf si{ln}lN) . . . z^ N{{ln}W \ (81) 

where Zk, (K = 1, 2, . . . , N), is the partition function of the /^-cluster, which obeys eqs.(31) 
and (32) or (33) (with N replaced by K), while A//r({/ n }|iV) denotes the {/ n }-realization 
dependent number of i^-clusters in an iV-chain with N nc non-catalytic sites. Evidently, for 
each realization {/„} these numbers J\fK({ln}\N) obey 

N 1 ({l n }\N) + 2N 2 ({l n }\N) + 3N 3 ({l n }\N)+ ... + N Af N ({l n }\N) = N (82) 

Note now that we have previously defined a -fT-cluster as being the set of all sites connected 
consecutively by arcs (see Fig.l). An equivalent definition of the X-cluster, which uses now 
the language of the intervals connecting consecutive non-catalytic sites is as follows: an 
interconnected K-cluster, whose partition function is determined by eqs.(31) and (32) or 
(33) (with iV replaced by K), is a subset of n (n < [(K — l)/2]) consequitive intervals 
l r +i, l r +2, lr+3, ■ ■ ■ , lr+n from the entire set {In}-, where all intervals a) are greater than unity, 
b) obey the conservation law / r+ i + / r+ 2 + . . - + l r +n = K — 1, and c) are necessarily bounded by 
two intervals (l r and l r+n+ \) of length unity. As we have already remarked, the latter condition 
implies that two pairs of non-catalytic sites appear at the extremities of the segment containing 
the K sites which automatically decomposes the chain into three independent parts. On the 
other hand, the condition that all intervals Ik, k = r + 1, r+2, . . . , r + n, are greater than unity 
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Figure 7: Example of a 10-cluster containing 4 inner intervals and 5 non-catalytic sites. 

insures that the -fT-cluster is interconnected and does not break up into smaller subunits (see 
Fig. 7). Consequently, we have from eq.(81) that 

N 

\nZ N ({l n }) = Mk({L}\N) hxZ K , (83) 

K=l 

and eq.(79) becomes 

N N 

(inz N (o) = p N - Nnc {i-p) Nnc Y. N ^ N ^\ N ) lnZ «i ( 84 ) 

N nc =0 K=l 

where now NK{N nc \N) reads 

N K (N nc \N) = J2^K({ln}\N), (85) 

and hence, defines the total number of K-clusters in all realizations of the A^-chain with a 
fixed number N nc of non-catalytic sites. 

The disorder-averaged pressure per site in the TV-chain with random, quenched placements 
of catalytic sites is then given by 

1 N 

(3P^en) {p) = £ Wr n{p) ln z ^ (g6) 

N^oo iv £ — * 
K=l 

where now wk,n(p) is the statistical weight of the AT-clusters in the TV-chain; wk,n{p) obeys 

N 

wk,n(p) = E P N ~ NnC ( J ~P) NnC N k{Nuc\N). (87) 

N nc =0 

Below we determine wk,n{p) explicitly. 

5.1 Calculation of the weights wk,n(p)- 

To fix the ideas, we start with the trivial case of (K = 1)- and (K = 2)-clusters. Consider a 
given realization {/„} of intervals. As may readily notice, a (K = l)-cluster, or "a free site" in 
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the terminology of the section 3, appears as soon as one has three consecutively placed non- 
catalytic sites. In other words, such a site appears as soon as any two consecutive intervals l r 
and l r+ \ are both equal to unity. Consequently, the number A/i({Z n }|-/V) of (K = l)-clusters 
in a given realization of the iV-chain with N nc non-catalytic sites can be written as follows: 

M({fn}W = X;<S(*r»l)*(fr+l,l) (88) 
r=l 

Then, using the definition of the lattice delta-function in eq.(10), we have that the total 
number N\{N nc \N) of {K = l)-clusters in all realizations of the iV-chain with a fixed number 
N nc of non-catalytic sites obeys: 

Nl (N nc \N) = EE^-^WH^E { ^ T^r^r-N+D = 

r=l {In} {In} JC T 

which yields, using the expansion 



(90) 



(i-r^-c-i n ^ 2 y Nnc _ 2 ) 

the following result: 

, s ( N - 2 \ f 1, if the N nc > 2, 

\ Jv nc — Z J [0, otherwise. 

Consequently, the weight wi j n(p) of the (K = l)-clusters is given by 

N nc =2 ^ nc lJ 

= N {l-pf + 2 (l-p)V (91) 



Next, we turn to calculation of W2,n describing the weight of the (K = 2)-clusters. Two such 
clusters, as we have already remarked, may only appear on the chain boundaries in the case 
when the sites i = 1 or i = N (or both) are catalytic, while two pairs of neighboring sites 
i = 2, 3 and i = N — 1,N — 2 are non-catalytic. Consequently, the number of (K = 2)-clusters 
in a given realization of an iV-chain with N nc non-catalytic sites is given by 

M 2 ({l n }\N) = S(h,2) S(l 2 ,l) + S(l Nnc+1 ,2) S(l Nnc ,l). (92) 

Hence, 

N 2 (N nc \N) = 1 E / - rE^ 1 ^ = -i / * fry^ t _ (n _ 2) 

I'nJ 
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and, making use of the expansion in eq.(90), we obtain 
N 2 (N nc \N) = 2 

Finally, we get 



N-3\ J 1, if the 2 < N nc < N — 1, 
N nc — 2/ 0, otherwise. 



w 2 , N = 2 (1 - pf p. (94) 

Now, in contrast to the very simple (K = 1)- and (K = 2)-clusters, clusters of larger size 
may be composed of several types of intervals. Let J\f^\{l n }\N) denote the number of K- 
clusters composed of n intervals in a given realization of an iV-chain containing exactly N nc 
non-catalytic sites. This number can be written down explicitly as 

^ ( k\{Iu}\N) = jfJ({l n }\K\N) + J$({l n }\K\N), (95) 

where 

n 

J$({l n }\K\N) = 2 (J[8(h >2)) S(l n+1 ,l)S(h + l 2 + ... +l n ,K) (96) 
i=i 

denotes the contribution from the if-clusters starting from either boundary site, i.e. "surface" 
.fT-clusters, while 

N nc -n n+r 

J$({l n }\K\N) = ]T ^fr' 1 ) ( II *(*i>2)) x 

r=l i=r+l 
X 6{l r+n+1 , 1) 8(l r+1 + l r+2 + ... + l r+ n + 1, K) (97) 

represents the contribution of the "bulk" if-clusters, i.e. such if-clusters which are entirely 
inside the chain and do not include any of the boundary sites. 

Summing M^\{l n }\N) over all the interval realizations obeying the conservation 
law in eq.(80), and next, performing in the result summation over all possible numbers n of 
sub intervals in a if -cluster, we find that for K / 1 and K / N, the total weight of if-clusters 
is given by 



p 



+ 



(98) 



wkAp) = P {K - 1)/2 (l-p) {K+3)/2 {2F K ^ ] l^ 

+ (l-p){N-K-l)F K - 2 (J^)}, 
while for K = N one has 

Details of these calculations are presented in Appendix B. 
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5.2 The thermodynamic limit N — > oo. 

Now, having calculated the weights wk,n(p) of if-clusters explicitly, we may rewrite eq.(86) 
as the sum of three contributions 

where the first term f3P[ quen \p) accounts for the contribution of (K = l)-clusters, 

PP^(P)= lim (-U^ln^Y (101) 

the second term denotes the contribution of a single spanning iV-cluster, 

(3P^ uen \p) = lim (±-w N , N lnZ N ), (102) 

while pp( quen ^ (p) takes into account all remaining possible if -clusters, 

pp ( g uen) {p) = fl g WK N{p) ln z \ {m 

~* °° V X=2 / 

In all these equations stands for the partition function of the corresponding if-cluster, 
which has been defined previously in eqs.(31) and (32). 

Now, using the results in eq.(91), we readily find that (5P[ quen \p) obeys 

j3P[ quen \p) = (1 -p) 3 ln(l + z). (104) 

Turning next to the contribution due to a single spanning iV-cluster, we have that it is given 
explicitly by 

x ^ ln C N + ln ( ^j== ) " N ln(r 2 )] ) . (105) 

Taking into account the explicit representation of the Fibonacci polynomials in eq.(35), we 
notice that for p > 1/2 their growth is suppressed by the exponentially vanishing factor 
(1 — p) N / 2 in the first line of eq.(105). On the other hand, for p < 1/2, the vanishing factor 

p N/2 in 

the first line of eq.(105) controls the large- A'" behavior. Consequently, as could be 
expected on physical grounds, the contribution of a single spanning iV-cluster is exactly equal 
to zero in the thermodynamic limit. 
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Consider now the form of (3P^ quen) (p) in eq.(103). Taking into account the explicit form of 
Zn in eqs.(31) and (32), and expanding Cn in eq.(32) in Taylor series in powers of (jiIt\) < 1, 
we may rewrite (5P^ quen \p) in eq.(103) as 

X=2 n=l v 17 / 

Evaluation of the sums entering eq.(106) is rather cumbersome and we present the details of 
such calculations in Appendix C. 

Taking into account both the contribution of the (K = l)-clusters in eq.(104) and that of 
in eq.(106) (see Appendix C) we arrive at the desired explicit expression describing 
the disorder-averaged pressure p(i uen ) (p) per site in the quenched disorder case: 



pP (quen \p) = (l-p) 3 ln(l + z)+p(5-7p + 3p 2 ) ln^ 1+ ^ + 1 ^ - p{l - P ^ ln(l + Az) + 

+ P {i- P fytlL ^lutl ( io7) 

which can be reformulated, (by expanding the denominator in elementary fractions), in the 
following form: 

0P(^)( P ) = (1 - pf ln(l + z) + p(5 - 7p + 3p 2 ) In ( V?T ^ + - P(1 ~ Vf m(l + 4s) - 
where the X± are given by 



p ± Vp(4 - 3p) . (109) 



2p(l - p) 

We note here that the first correction to the thermodynamic limit result in eqs.(107) or (108), 
(which shows how fast the thermodynamic limit is approached with respect to the chain's 
length), should be proportional to the first inverse power of N, as follows from the expansion 
eq.(C.12) (see Appendix C). 

We now consider the asymptotical behavior of p(i uen )(p) for different limiting cases. In 
the asymptotic limit p^Owe find from eq.(107) that 

0P {quen) (p) = In (1 + z) + In ( 1 ^ 3Z+ J 2 ) P + 0(p 2 ), (110) 

\ (1+2) / 
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0.2 0.4 0.6 0.8 1 
P 

Figure 8: Mean particle density n versus mean density p of catalytic sites for different values of 
z = exp(/3/x). Curves with signs (ann) and (quen) depict the behavior of the mean density for annealed 
and quenched random distributions, respectively. 



where the first term represents the Langmuir pressure, and the second one - the small-p 
correction to it. Note that already the first correction term differs significantly from the first 
correction term to p( ann \p) found in the annealed disorder case, eq.(61). Next, in the limit 
p ~ 1, we find that p(i uen ) obeys 

/,/tti; i i\ I (i + Vi + 4z N ) 2 \ 

ppiiuen) {p) = ln +^+ j + ln V 4 ^__ ; (1 - pf + O ((1 - pf) , (111) 

in which expansion the second term is also different from the one obtained in the annealed 
disorder case, eq.(63). 

Now, turning to the analysis of the large- z behavior of p(i uen ) we notice that the behavior 
differs completely for p = 1 and for p < 1, which signifies that here, (as in the annealed 
disorder case), p = 1 is a special point. In the case p = 1 we have that on the righthand side of 
eq.(107) all terms except for the second one vanish and hence, p(i uen \p) = p( re 9)(L = lor 2), 
eq.(36). Consequently, for p = 1 the large-z behavior of P^i uen )(p) obeys eq.(76), similar to 
p{re 9 )^ L = i or2 ) and to P^ ann \p), which is, of course, not surprising. On the other hand, 
for p < 1 the large-z behavior is rather different from the behavior observed in the annealed 
disorder case. Here, we find for z > (1 - p)~ 2 that p(<?« en )(p) obeys 

PP {quen) {p) = ^ + ^ ln (z) +0(1), (112) 
i.e. in the quenched disorder case the prefactor in the leading z-term depends on p, while in 
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Figure 9: Mean particle density n in the annealed and quenched disorder cases, eqs.(66) and (113), 
versus the activity z for different values of p. The notations are the same as in Fig. 8. 



the annealed disorder case this prefactor was found to be independent of p, which caused a 
rather strange behavior of the mean particle density. 

Differentiating eq.(108), we find that the particle density in the quenched disorder case is 
explicitly given by 

z ( 2 \ 2z 2p(l-p) 2 z 



n (9«en)(p) = (l_p)3_A_ + ^5_7 p + 3p 2^_ 



+ 4z + \/TTIz 1 + 4z 

72 \ m+4 



V p(4-3p)(l + 4,) (l + v TT4l) 2 ^oU? X-) ^_^_^y + J 

We note that again this result differs considerably from the mean-field prediction n = (eq. (8) 
with K. = oo). 

From eq.(113) we find then that the asymptotic behavior of rS- quen \p) in the small-p limit 
obeys 

„(^ M = ^- (i+ ^;;) +z2) P+0( ^, (114) 

while in the limit p ~ 1 it follows 

(i uen \p) = n( re9 \L = lor2) + - . ^ . (1 - pf + O ((1 - pf) , (115) 

y > V ' {l + Az) (l + 2z + Vl + 4z) V ! 



n 



in which equations the first corrections to the Langmuir and the regular cases depend very 
differently on the activity z when compared to the annealed disorder case. Note also that, 
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as depicted in Fig.8, for any fixed z the mean particle densities in the annealed and in the 
quenched disorder cases show a completely different behavior as functions of the mean density 
p. The difference becomes more pronounced with increasing z. 

Note that the catalytic efficiency in the annealed disorder case turns out to be lower than 
in the quenched case, as can be inferred from the fact that the A particle density is always 
higher in the former case (see Fig. 9). 
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Figure 10: Compressibility kx as a function ofp for the annealed and for the quenched disorder cases. 



Finally, we find that as z — > oo, the mean particle density n^ quen ^ tends to 

lim n^ uen \p) = 1 V —, (116) 

z->oo 1 + p 2 

which contradicts apparently the behavior observed in the annealed disorder case, where we 
found that \\T0. z ^ oo n^ ann \p) = 1, regardless of the value ofp, (provided that p < 1); it also 
differs from our predictions for the case of a regular placement of catalytic sites, for which 
lirn^oo n( re f)(L) = 1 - p for p < 1/2 and lim^oo = 1/2 for p = 1/2 and p = 1. Note 

also (despite of the fact that here p = 1 also appears as a special point in regard to the large-z 
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behavior of the disorder-averaged pressure) that here, contrary to the annealed disorder case, 
n( guen \p) does not show any discontinuity in the limit p — ► 1 at z = oo. 

To close this final section we discuss the behavior of the compressibility in the quenched 
disorder case. From eqs.(108) and (113) we find that in the small-p limit the compressibility 
follows 

z (1 + 3z + z z ) 

while in the opposite limit p ~ 1 it is described by a more complicated expression of the form: 

r '4-' = = 1 or 2) + ^g||g^ (1 - P f + O (d - P)*)H8) 

It follows then that the coefficient in the term proportional to (1 — p) 3 is positive only for 
z < 2 and negative for z > 2. This implies, in view of the discussion presented at the end of 
the previous section, that for z < 2 the compressibility k^ uen ^ is a non-monotonic function 
of p. On the other hand, for z > 2 the compressibility k^ uen ^ seems to be always increasing 
with p. This is distinct from the behavior found in the annealed disorder case when k^ nn ^ is 
non-monotonic function for any z. 

In Figs. 10 we depict the behavior of the compressibility as a function of p in the quenched 
disorder case and compare it to the behavior observed in the annealed disorder case. These 
figures suggest that for z < 2 the compressibility k^ uen ^ is a non-monotonic function of p, 
while for z > 2, (as exemplified here by the case z = 5) it is monotonic. We also find that 
the most dramatic difference between the three different modes of placing the catalytic sites 
is seen in kr- 



6 Conclusion. 

In this paper we have studied the properties of the catalytically-activated annihilation A+A — ► 
reaction on a one-dimensional lattice, in which some lattice sites possesses special "catalytic" 
properties; reaction takes place when at least one of two neighboring adsorbed A particles, 
undergoing continuous exchanges with a particles reservoir, sits on a catalytic site. 

We have focused here on three different situations. First, we have considered the case 
when the catalytic sites are placed periodically, forming a regular sublattice, in which case we 
obtained the exact solution in a straightforward manner. Next, we turned to the disordered 
case and studied the reaction properties for both annealed and quenched randomness in 
the distribution of the catalytic sites. We have shown that in the annealed disorder case 
the model reduces to a one-dimensional lattice gas with an effective three-particle repulsive 
interaction. We have developed a combinatorial formulation of the model which allowed us to 
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obtain the exact solution. Next, we have demonstrated that in the (most complex) situation 
with quenched disorder the problem of computation of the average logarithm of the partition 
function can be reduced to the problem of the enumeration of all possible interconnected 
clusters in finite lattices with a fixed number of catalytic sites. We have calculated the 
weights of these clusters exactly, using a combinatorial procedure, and we have found an 
exact expression for the disorder-averaged pressure. 

Apart from the results on the disorder-averaged pressure, we have determined exact 
asymptotic expressions for the mean-particle density and for the compressibility. We have 
shown that the behavior of these properties is substantially different in systems with an- 
nealed and with quenched disorder. Both differ considerably from the mean-field result of 
eq.(8). Furthermore, we have observed that in systems with annealed disorder the mean 
particle density tends to unity in the limit when the chemical potential fi tends to infinity, 
and this for any limited catalytic sites density. On the other hand, we have established that 
the mean particle density in the quenched disorder case tends to a finite value, namely to 
(1— p + p 2 )/{l+ p 2 ) < 1 when u tends to infinity. As well, we have demonstrated that in 
the annealed disorder case the compressibility appears to be a non-monotonic function of p 
for any /i, while in the quenched disorder case the compressibility shows a non-monotonic 
behavior as a function of p only for /i < n C rit = P 1 m (2). 

We close by noting that the model studied here furnishes another example (see, e.g., 
Refs. [19, 36-38]) of a ID Ising-type system with random multisite interactions which admits 
an exact solution. 
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Appendix A 



The partition function in eq.(55) can be represented as 

(Z N (0) = h - h, (A.l) 



where 



/ N N + l 

2m J c t \ p ) 

- + z{l - P))r ~ r 2 + t*Y\ (A.2) 



and 



x 



2mp Jc >\ (l-p)(l- r ) ; 

- <Iz£> (i + z(i _ P ))r - r 2 + r 3 ) \ (A.3) 
V p p J 



To evaluate the integrals in eqs.(A.2) and (A.3), let us first express the denominator of 
the integrands in terms of elementary fractions; this gives 

1 1 



^-P? _ (1~P) (1 + z{l _ p))T _ T 2 + T 3^ " ~ ' i (/:! /->)(/ : - /, )( i - / r T r" 
1 



hfa - ti)(t 2 - *s)(l - tjV) t 3 (t 3 - i 2 )(ia - «i)(l - ii" V) ' 
where ti, i 2 and i 3 are three roots of the cubic equation 



(A.4) 



t 3 _ t 2 _ 0_£l (i + 3(1 _ p))/ + z(1 y) = 0. (A.5) 
p p 

Expanding next (1 — t~ l r) in a Taylor series with respect to r and taking advantage of the 
definition of the lattice delta-function in eq.(10), we find that I\ of eqs.(A.2) is given by 

j.2 



( ptj-ph-il-p) ( z(l-p) \N +1 

l p(h - t 2 )(t! - 1 3 ) v ti / 



P^-pt 2 -(l-p) p(l-p) \^v+i 

p(t 2 - h)(t 2 - 1 3 ) v t 2 / 

P* 3 - - (1 - P) /^(l - P) ^ JV+1 



p(t 3 -t2){t 3 -h) V t 3 

On the other hand, the function in eq.(A.4) and also the expression 

1 -p-pT+(l -p)T 2 \N+2 



}• (A.6) 



( i-p-pT + {i-p)T- yy+^ 

y (l-p)(l-r) ) (A - 7) 
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are both analytic functions of the variable r. Hence, in virtue of eq.(10), it follows that the 
integral in eq.(A.3) equals zero and hence, I2 = 0. 

Furthermore, in terms of the auxiliary functions R and Q, Eqs.(57) and (58), the roots of 
the cubic equation (A. 5) can be written as follows [46]: 



1 

3 + 



(r + VQ 3 + R 2 ) 1/3 + {r- VQ 3 + R 2 ) 



1/3 



(A.8) 



t 2 



l l r 

3 ~ 2 



1/3 



R + VQ 3 + R 2 ) +[R- VQ 3 + R 2 

1 /3 

- (r - VQ 3 + 



1/3- 



1/3 



(A.9) 



and 



+ 



1 1 

3 ~ 2 



(r + VQ 3 + R 2 ) 1/3 + - VQ 3 + R 2 ) 1/3 

i? + vq 3 +^ 2 ) 1/3 -(r- vq 3 + ^ 1 



+ 



(A.10) 



Note next that the characteristic sum + R 2 is less or equal to zero for any value of the 
parameters p and z; hence, all three roots of the cubic equation (A. 5) are real. Noticing also 
that the ratio R/\/—Q 3 is bounded, — 1 < R/^/—Q 3 < 1, we find that the roots can be 
expressed in a more convenient fashion as: 

R 

4- 2\ - -( ) cos ! — ni'ccos 1 — 



h = ^ + 2,/^Q cos (J arccos ( ^ r -^- 3 

1 — T /I / /">' 

t 2 = 3-2v^ S m(- 



■ arcsm 



and 



(A.ll) 
(A.12) 



t 3 =\- 2^Qsm g + I arccos (^=)) 



Noticing now that all roots tj, (i = 1, 2, 3), of eq.(A.5) obey 



pt 2 _ pt . _(!_„) = 2(1 - p) - 1) (^^) , 



while 



t\ — £2 = 2 a/— 3Q cos f — + - arccos 
V 6 3 



2t,; 
i? 



ti — ts = 2^/-3Q cos — - arccos ( 



and 



*2 — *3 = 2\/— 3Q sin f- arccos 
V 3 



v 3 ^ 



(A.13) 

(A-14) 

(A.15) 
(A.16) 

. , (A.17) 

V^Q 3JJ 

and substituting the results in eqs.(A.ll) to (A.13) into eq.(A.6), we find, eventually, the 
result in eq.(56). 
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Appendix B 

Summing ({l n }\N) over all the interval realizations {l n } obeying the conservation law in 
eq.(80), and using the integral representation of the Kronecker function in eq.(10), we obtain, 
after some straightforward calculations, that: 



{'»} 

while 



=<««-> ( K „"-r) (%----*)■ < B - 2 » 

{In} 

Performing next the summation over all possible numbers n of subintervals in a -fT-cluster, 
we find that the total number Nx(N nc \N) of if -clusters for all possible realizations of an 
TV-chain containing a fixed number N nc of non-catalytic sites is given by 

[K/2] fK-l-n\ fN-K-l\ 

+ 



+ 5(K,N)[2Y J S(n,N nc )( ^ _ ^ n J + 

n=l ^ 7 

[(^-l)/2] ,„ „ v [{K+l)/2] , . 

+ j: s(n, Nnc - 1) r j - n ) + e + 1) r) } 

n=l n=l ^ 7 



N - 2 

+ <5(if,l) (1 - <5(iV nc , 1)) (l-<5(iV nc ,0)) iV nc ( ). (B.3) 



N nc -2 



Now, representing the weights wk,n{p) of JT-clusters as 

wk,n (p) = (p) + u-g^ (p) , (B.4) 
where w KN (p) and w KN (p) denote the weights of the "bulk" and "surface" if-clusters, 



respectively, we find, summing over all possible numbers N nc of non catalytic sites, that these 
weights are defined explicitly by 

[(K-l)/2\ /Tjr n N oo /-, x JV n 

«■<,..,- 

n=l x 7 Ar nc =0 

and 



n=l V 7 Af nc =0 V ^ 7 V 7 
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Noticing next that 

AT, 



AT nc =0 

we obtain that w^ N (p) obeys 

^(P) = (1-P) 2 !^- 1 £ [(l"P)(iV-K-l) + 2]( ^^p) (^) , (B.8) 

and consequently, in virtue of the definition of the Fibonacci polynomials, eq.(35), the total 
weight of the " bulk" If -clusters can be expressed by 

W W N (p) = p (K-l)/2 {1 _ p) (K + 3)/2 x 

x [(l- p )(N-K-l) + 2]F K - 2 (^^. (B.9) 
In similar fashion, we find that the total weight of "surface" if-clusters w^ N (p) obeys 

«,<*>) = 2 P ^ (1 - P )^^ F K . X [s^) • (B.10) 
Combining these results, we arrive eventually at the general formulae in eqs.(98) and (99). 

Appendix C 

In order to evaluate the limiting behavior of the rather complicated sums entering eq.(106), 
it is expedient to introduce an auxiliary generating function of the form: 

oo 

F{M = Y,Sn(Ot n , (C.l) 

N=3 

where 

M0 = £^^. (C2) 

K=2 

Once JF(£, t) and Gn{Q are determined, one obtains (3P^ quen \p) directly. As one may verify 
readily, pP^ip) 

in eq.(106) can be expressed in terms of Gn(0 as 



lnr 2 



n=l 
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We turn now to the calculation of the generating function in eq.(C.l). Substituting the explicit 
form of wk,n(p), eq.(98), into eq.(C.2), and this in eq.(C.l), and interchanging then in the 
final result the order of summations (over K and N), we find that r) can be represented 
as a sum of two components, 

F(S,t) = F 1 (S,t)+F 2 (S,t), (C.4) 

where 

y K=2 vv ■ K/ N=K+1 

(C.5) 

and 



, (1~P) 3/2 

' pl/2 

Using next the evident integral equality 



oo TV 



1 

TV 



1 

dvw"- 1 , (C.7) 



as well as the explicit expression of the generating function of the Fibonacci polynomials [40] 

oo 

we obtain the following integral representations of iq(£,r) and i*2(£,r): 

Fi«,r) = (i-p)V( 3 Ttt^t " V ( c - 9 ) 



and 



^«,r) = 2(1-„)Y'V€ / 1 n % 7 Y - (CIO) 

■A) (1-rv) (l - port - p(l - p)(vt£)*) 

Now, the asymptotic behavior of behavior of Gn{£) as N — ► oo can be deduced, in a standard 
fashion, by analysing the critical behavior of the generating function r) in the vicinity 
of the singularity closest to the origin [47]; that is, here, in the vicinity of r = 1. Expanding 
r) in the vicinity of the singular point r = 1, we obtain 

^^' rj i_ K -p(i-p)e 2 l-r + 

2(1 - n)V /2 r 3 £(l -p + f ) 
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Next, in virtue of the Tauberian theorems [47], it follows that Gn{Q exhibits the following 
asymptotical behavior as N — ► oo: 



Gn(0 



p(l-p) 4 e _ ( 2(l-p)V/ 2 r 3 e(l-p + g) ] 1 /l, ( . 
1 -p^-p(l -p)? 2 I I - p£ - p{l - p)t 2 IN °\N.' ' 



and consequently, we find from eq.(C.12) that the particular values of the function Gn(£.) 
entering eq.(C.3) are given explicitly by: 



&o(£ = i) =p(i-p) ; 



2. <%U0 



9£ 



p (p 2 -3p + 3), 



and 



5oo U 



T2 



p(l - p) 4 (-r 2 /ri) 



3« 



(C.13) 



(C.14) 



Tiy / l_ p( _l)n (r2/rir _ p(1 _ p)(r2/ri) 2n- 

Substituting the results in eqs.(C.13) and (C.14) into eq.(C.3), we thus find that f3P^ quen \p) 
obeys 

ppt en \ P ) = p(5 - 7p + 3p 2 ) In ( ^f +1 ) - 5(lzrt! ln(1 + 4z) + 



+ p (i-p)4yLill 

jfe « l-M-l^^/rir-Ml-^^/n) 2 " 
which, in combination with the expression for f3P[ quen \p), eq.(104), leads to eq.(107) 



5n 



(C.15) 
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